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Abstract 

We study the distribution of the exponential functional /(^,??) = exp{^t-)dr]t, where 
^ and T] are independent Levy processes. In the general setting, using the theory of Markov 
processes and Schwartz distributions, we prove that the law of this exponential functional sat- 
isfies an integral equation, which generalizes Proposition 2.1 in [9]. In the special case when 
r/ is a Brownian motion with drift, we show that this integral equation leads to an important 
functional equation for the Mellin transform of 1(^,77), which proves to be a very useful tool 
for studying the distributional properties of this random variable. For general Levy process ^ 
(rj being Brownian motion with drift) we prove that the exponential functional has a smooth 
density on M \ {0}, but surprisingly the second derivative at zero may fail to exist. Under the 
additional assumption that has some positive exponential moments we establish an asymp- 
totic behaviour of P(/(^,?7) > x) as a; — t- -|-oo, and under similar assumptions on the negative 
exponential moments of ^ we obtain a precise asymptotic expansion of the density of /(^, r/) as 
x — >■ 0. Under further assumptions on the Levy process ^ one is able to prove much stronger 
results about the density of the exponential functional and we illustrate some of the ideas and 
techniques for the case when ^ has hyper-exponential jumps. 
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1 Introduction 



In this paper, we are interested in studying distributional properties of the random variable 



where ^ and 77 are independent real-valued Levy processes such that ^ drifts to —00 and IE[|'Ci|] < 00 
and IE[|r7i|] < 00. 

The exponential functionals /(^, rj) appear in various aspects of probability theory. They describe 
the stationary measure of generalized Ornstein-Uhlenbeck processes and the entrance law of positive 
self-similar Markov processes, see [6, 9]. They also play a role in the theory of fragmentation processes 
and branching processes, see [4, 22]. Besides their theoretical value, the exponential functionals are 
very important objects in Mathematical Finance and Insurance Mathematics. They are related to 
Asian options, present values of certain perpetuities, etc., see [10, 17, 14] for some particular examples 
and results. 

In general, the distribution of exponential functionals is difficult to study. It is known explicitly 
only in some very special cases, see [8, 14, 19]. Properties of the distribution of I{^,rj) are also of 
particular interest. Lindner and Sato [26] show that the density of /(^,?7) doesn't always exist, and 
in the special case when ^ and 77 are specific compound Poisson processes, distributional properties 
of /(^,?7) can be related to the problem of absolute continuity of the distribution of Bernoulli con- 
volutions, which dates back to Erdos, see [12]. The distribution of 1(^,77), when C,s = —s and in 
some other instances, is known to be self-decomposable and hence absolutely continuous, see [5, 18]. 
When ?7 is a subordinator with a strictly positive drift, the law of the exponential functional I{C,,'r]) 
is absolutely continuous, see Theorem 3.9 in Bertoin et al. [5]. Some further results are obtained in 



The asymptotic behaviour F{I{C,,r]) > x), as x — 00, is a question which has attracted the 
attention of many researchers. In the general case, but under rather stringent requirements on the 
existence of exponential moments for ^ and absolute moments for 77, it has been studied in [25]. The 
special case when rjt = t has been considered in [27, 31, 32] and properties of the density of the law 
of I{^,ri) at zero and infinity have been studied by [19, 21, 28] and results such as asymptotic and 
convergent series expansions for the density have been obtained. 

The first objective of this paper is to develop a general integral equation for the law of I{^,r]) 
under the assumptions that IE[|^i|] < 00, ]E[^i] < 0, IE[|?7i|] < 00 and ^ being independent of 77. 
Using the fact that in general /(^, 77) is a stationary law of a generalized Ornstein-Uhlenbeck process, 
Carmona et al. [9] show that if has jumps of bounded variation and rjt = t then the law of 
/(.^, rf) satisfies a certain integral equation. We refine and strengthen their approach and using both 
stationarity properties of /(^,?7) and Schwartz theory of distributions, we show that in the general 
setting the law of /(^,?7) satisfies a certain integral equation. This equation is important on its own 
right, as demonstrated by Corollary 1, but it is also amenable to different useful transformations as 
can be seen from the discussion below. 

The second main objective of the paper is to study some properties of I^^^ := rj) in the specific 
case when r]s = fis + aBg, where Bg is a standard Brownian motion. Quantities of this type have 
already appeared in the literature, see [14], but have not been thoroughly studied. The latter, as it 
seems to us, is due to the lack of suitable techniques, which are available in the case when 77^ = s, 
and in particular due to the lack of any information about the Mellin transform of /(^, rj), which is 
the key tool for studying the properties of I{^,r]), see [19, 21, 27]. We use the integral equation (2.3) 




(1.1) 



[24, 29, 30, 35]. 
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and combine techniques from special functions, complex analysis and probability theory to study the 
Mellin transform of J^^o-; which is defined as Ai{s) = E [(/^ ^>o}] • In particular we derive 
an important functional equation for Ai{s), see (3.13), and study the decay of Ai{s) as Im(s) — oo. 
These results supply us with quite powerful tools for studying the properties of the density of I^^a- via 
the Mellin inversion. Furthermore, the functional equation (3.13) allows for a meromorphic extension 
of Ai{s) when C, has some exponential moments. This culminates in very precise asymptotic results 
for P (/^,o- > x), as X — >■ oo, see Theorem 5, and asymptotic expansions for k{x), the density of 
J^^o-; as X — )■ 0, see Theorem 4. The latter results show us that while k{x) G C°°(]R \ {0}), rather 
unexpectedly k"{0) may not exist. Finally, we would like to point out that while the behaviour of 
P > x), as X — )■ oo, might be partially studied via the fact that J^^o- solves a random recurrence 
equation, see for example [25], the behaviour of k{x), as x — >■ 0, seems for the moment to be only 
tractable via our approach based on the Mellin transform. 

As another illustration of possible applications of our general results, we study the density of J^^g- 
when ^ has hyper-exponential jumps (see [7, 8, 20]). This class of processes is quite important for 
applications in Mathematical Finance and Insurance Mathematics, and it is particularly well suited 
for investigation using our methods due to the rich analytical structure enjoyed by these processes. 
In this case we show how to derive complete asymptotic expansions of k{x) both at zero and infinity. 
We point out that our methodology is not restricted to this particular case, and can be easily applied 
to more general classes of Levy processes. 

The paper is organized as follows: in Section 2, we study the law of rj) for general independent 
Levy processes ^ and r] and derive an integral equation for the law of 77); in Section 3, we specialize 
the results obtained in Section 2 to the case when r]s = ns + aBg and, employing additionally various 
techniques from special functions and complex analysis, we study the properties of the density of I^^a- 
Section 4 is devoted to some applications of the results derived in the previous section. In particular, 
we study the asjTiiptotic behaviour at infinity of the tail of /^^o- and of its density at zero, and in 
the case of processes with hyper-exponential jumps, we show how these results can be considerably 
strengthened. 



2 Integral equation satisfied by the law of /(^, 77) 

Let us introduce some notation which will be used throughout this paper. The main underlying 
objects are two independent Levy processes ^ and 77 defined on a probability space (^2, J-", P). As is 
standard, we assume that both processes are started from zero under the probability measure P. 

Assumption 1. Everywhere in this paper we will assume that 

E[|ei|]<oo, E[ei]<0, E[|77i|]<oo. (2.1) 

The characteristics of the Levy processes ^ and 77 will be denoted by (6^, cr^, 11^) and (6^, a^, 11^). 
In particular Tl^^idx) and 11^ (dx) are the Levy measures of ^ and 77, respectively. We use the following 
notation for the double-integrated tail 

=(+) f°° =(-) f°° 

(x) = / Ii^{{y,oo))dy and (x) = U^{{-oo, -y))dy, 

J X J X 

=(+) =(-) 

and similarly for 11^ and 11^ . Using the Levy-Ito decomposition (see Theorem 2.1 in [23]) it is 
easy to check that Assumption 1 implies that the above quantities are finite for all x > 0. 
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We define the Laplace exponents ip^{z) = In (E [e^^^]) and 'ipniz) = In (E [e^''^]), where without 
any further assumptions and ipj^ are defined at least for Re{z) = 0, see [3, Chapter I]. The Laplace 
exponent can be expressed in the following two equivalent ways 

^jj^(z) = -^z^ + b^z+ / {e''' -1- zx)U^{dx) (2.2) 

1-2 / poo I I N /"oo 



-^z^ + b^z + zUj {w)e^'dx + J {x)e-'''dx\ 



with a similar expression for iprj- The first equality in (2.2) is essentially the Levy-Khintchine formula 
(see Theorem 1 in [3]) with the cutoff function h{x) = 1. The standard choice for the cutoff function 
in the Levy-Khintchine formula would be l{|a;|<i}, however it is well-known that if E[|,^i|] < oo then 
we can take a simpler cutoff function h{x) = 1. The second equality in (2.2) follows easily by 
repeated integration by parts. Note that according to (2.2), we have = ip'^{0) = E,[^i] and similarly 

6, = V^;(o) = EM. 

We recall that the exponential functional I{^,r]) is defined by (1.1), its law will be denoted by 
■m{dx) := P(/(^,?7) G dx). The density of /(^,?7), provided it exists, will be denoted by k{x). 

Our main result in this section is the derivation of an integral equation for the law of I{C,,ri). 
This equation will be very useful later, when we'll derive the functional equation (3.13) for the Mellin 
transform of the exponential functional in the special case when 77 is a Brownian motion with drift. 
The main idea of this Theorem comes from Proposition 2.1 in [9]. 

Theorem 1. Assume that condition (2.1) is satisfied. Then the exponential functional I {^,7]) is well 
defined and its law satisfies the following integral equation: for f > 

j m(dx)^ dv + -^vm{dv) + {^j ^ -^^ {^^ ~) m(dx)^ dv + ^In — ^ m(dx)^ dv 

+ {- I n^^V - x)m{dx)] dv + (- I u[\x - v)m{dx) ] dv (2.3) 



/•oo ^ f^=(+) \ / 1 f°°={-) \ 

J — J {w — x)m{dx)dw j dv — I J ~ J ~ 'w)m{dx)dw j dv = 0, 

where all quantities in (2.3) are a.e. finite. Equation (2.3) for the law of —rj) on (0, 00) describes 
m(dx) on (—00, 0). 

The proof of Theorem 1 is based on the so-called generalized Ornstein-Uhlenbeck (GOU) process, 
which is defined as 



Ut = Ut{^,ri) = xe^* +e^' f e ^'~dr]s = xe^* + f e^=-d?7^, for t > 0. 

Jo Jo 



(2.4) 



Note that the GOU process is a strong Markov process, see [9, Appendix 1]. Lindner and Mailer 
[25] have shown that the existence of a stationary distribution for the GOU process is closely related 
to the a.s. convergence of the stochastic integral e^^"d?7s, as t — )■ 00. Necessary and sufficient 
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conditions for the convergence of I{^,f]) were obtained by Erickson and Mailer [13]. More precisely, 
they showed that this happens if and only if 



lim C,t = —oo a.s. and 



t—>oo 



-e,e] 



i + /;°^'''^'n5(M\(-^,^))d^ 



n^(dy) < oo. (2.5) 



It is easy to see that Assumption 1 implies (2.5). Hence is well-defined and the stationary 

distribution satisfies Uoo = I{C,,f])- This identity in distribution is the starting point of the proof of 
Theorem 1. 

As the proof of Theorem 1 is rather long and technical, we will divide it into several steps. We 
first compute the generator of U, here denoted by This result may be of independent interest, 
therefore we present it in Proposition 1 below. Then we note that the stationary measure m(dx) 
satisfies the equation 

poo 

il^^'^f{x)m{dx) = 0, (2.6) 



where / is any infinitely differentiable function with a compact support in (0,oo). Indeed, (2.6) 
follows from (2.1) in [9] or from the definition of infinitesimal generator and the observation that, for 
all t > 0, 

poo POO 

E[f{Ut)]m{dx) = / f{x)m{dx). 
Jo 

Finally, an application of Schwartz theory of distributions after rephrasing (2.6) gives (2.3). 

We start by working out how the infinitesimal generator of U, i.e. H*-^-*, acts on functions in 
/C C Co(M), where 

/C = {fix) : fix) e C,2(M), fie^ G CHR) n Co(M)} 

n { fix) = 0, for X < 0; /'(O) = /" (0) = 0} (2.7) 

and C^(M) stands for two times differentiable, bounded functions with bounded derivatives on M and 
Co(M) is the set of continuous functions vanishing at ±oo. Denote by and ( resp. and 
D^) the infinitesimal generators (resp. domains) of ^ and t]. Note that 

^^^^fix) = kf'ix) + ^f"ix)+ [ ifix + y)-fix)-yrix))U^idy) (2.8) 



= hf'ix) + ^f"ix)+ f"ix + w)U^ iw)dw+ f"ix-w)U^ iw)dw, 

with a similar expression for The first formula in (2.8) is a trivial modification of the form of the 
generator of Levy processes for the case when the cutoff function is /i(x) = 1, see [3, p. 24], whereas 
the second expression follows easily by integration by parts, the fact that / G /C and lE[|^i|] < oo. 
Finally, we are ready to state our result, which should strictly be seen as an extension of Proposition 
5.8 in [9] where the generator 2.^^^ has been derived under very stringent conditions. 
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Proposition 1. Assume that condition (2.1) is satisfied. Let / G /C, g{x) := and (j){x) :- 

/(e^')- Then, / G 2)'', G D« and 



,2 



a 



k9{x) + -4-^9'' 



X] 



g'{v)T{\ ^ In ^) dt; + j g\v)T]:^^ ( In ^) 



cr, 



.=(+) 



+ ^r(x) + / r(x + «;)n^ Hdw;+ / r(x-ti;)n^ Hd^. (2.9) 
^ Jo Jo 

Proof. The main idea is to use the definition of the infinitesimal generator and Ito's formula. Let 
/ G and note that by definition 



= Urn "-'^'^'"-^W^ li.nl IE 



Using the fact that (f/t)t>o ^ semimartingale and / G /C, we apply Ito's formula to f{Ut) to obtain 

f{U,)-f{x)= f f{U,.)dU, + \ f f"{Us^)d[U,Ur, + Y,{f{Us) - f{Us^) - AUsfiUs^)) . (2.10) 
Jo ^ Jo 

Now, let Ht := e^* and Vt := x + e~^^'dr]s, and note that f/t = HtVt. Hence by integration by parts 

Ut = x+ [ Hs^dVt+ [ Vs^dHs + [H,V]t. 
Jo Jo 

Using the Levy-Ito decomposition (see Theorem 2.1 in [23]) and Assumption 1, we find that the 
Levy processes ^ and r] can be written as follows 



= (T^Bt + b^t + Xt, Vt = (^r,Wt + \t + Yu 



(2.11) 



where B and W are Brownian motions, X and Y are pure jump zero mean martingales, and the 
processes 5, W,X and Y are mutually independent. Then we get 



Vt = x 



br, [ e-^'-ds + a^ [ e'^'-dlVs + Nt, 
Jo Jo 



where A^^^ = e ^""dy^ is a pure jump local martingale. On the other hand using Ito's formula, we 
have 

Ht = e^^ = l+ f e^'-dis + \ f e^'-d% i]l + ^ e«=- (e^«^ - A^, - 1) 
Jo ^ Jo 

= l+(b^ + ^\ I e^^-ds + f e^^-dB, + iV^ + J] (e^«= - A^, - l) , 
\ /Jo Jo 

where A^, = e^^dX^ is a pure jump local martingale. Therefore, we conclude that 



I e 





-«=-ds„(T^ / e-«^-diy. 



+ ^AVsAHs = a.s., 

* s<t 



since AVg = e ^"-Arjs, AH^ = ifs_(e^^= — 1) and the fact that ^ and r] are independent and do not 
jump simultaneously a.s. This implies that 

Ut = x+ [ Hs-dVs+ [ Vs-dHs = x+ [ e^^-dVs+ [ Vs-dHs. 
Jo Jo Jo Jo 

Using the expressions of H and V, we deduce that 



Ut = x + b„t + arjWt + I e^'-dN, + { + 



V.^e^'-ds 



X 



Jo Jo ^ ^ 

+ Kt + + b,t + (^b^ + Us-ds + J2 Us- (e^«» - - l) , 



where 



e^'-dN, 



K-diV, = Yt 



= ar^Wt + a^ [ U^.dB,. 
Jo 

From the definition of K and K"^, and the mutual independence of B, W, N and A^, we get for the 
continuous part of the quadratic variation of U 

[U, UYt = \K\ K% = alt + al f U^ds. 

Jo 

Putting all the pieces together in identity (2.10), we have 

fiUt) - fix) = Mt + b, f'iUs-)ds + (^b^ + y) I f'(Us-)Us-ds 

+ E f'iUs-)Us- (e^«» - Ae. - l) + ^ £ f"{Us-)ds + nUs-)Ulds 



s<t 



where M is a local martingale starting from and M describes the integration with respect to K 
and K'^ in the expressions above. Using the fact that f & K, implies f{x) = for x < and 
x\f'{x) \ + x^\f"{x)\ < C{f) < oo, we deduce that Mt is a proper martingale as all other terms in 
the expression above have a finite absolute first moment. Furthermore applying the compensation 
formula to the jump part of f{Ut) we get 



E 



s<t 



E 



[ nUs.)Us(l - y - l)U^{dy)] ds 
Jo \Jym J 
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Similarly, using the fact that AUs = Arjs when Arjg ^ and AUg = Us^^e^^" — 1) when A.^^ 7^ 
(see the definition of U) we get 



E 



^(/(f/.)-/(f/._)-A[/,_/'(f/,_)) 



s<t 



E 



if (Us- + z)- /(f/,_) - zfiUs-)) U,{dz)ds 



Jz£ 



+E 



Jo Jy£R ^ ^ ' ^ . 



Finally, as / G /C, we derive 



E 



/(f/Oj - /(x) = 6,E 
2 
2 

E 



^E 



r(f/._)f/f_d. 



65 + ^lE 



/'(t/,_)f/,_ds 



L^'O 



2 



r\Vs-^)ds 



VJO 



VJo 

t 



E 



r / (/(f/,_ + ^) - /([/,_) - z/'(?7,_))n,(d^)ds 

r / (/(t/.-e^) - /(f/.-) - i//'(f/._)f/._)n5(d|/)d. . 
and dividing by t, letting t go to and recalling that Uq = x a.s., we obtain for / G /C the identity 

(2 \ 2 2 

+ / ff{x + z)-f{x)-zf'{x))u,{dz)+ [ (fixey) - fix) -yxf'{x))u^{dy),{2.12) 



and therefore the infinitesimal generator of U satisfies 

In order to finish the proof one only has to apply integration by parts. 
The following Lemma will also be needed for our proof of Theorem 1. 



□ 



Lemma 1. Assume that condition (2.1) is satisfied. Let u{dv) denote the measure in the left-hand 
side of formula (2.3). Then |z/|(df) and hence u{dv) define finite measures on any compact subset of 



(0, 00) and for any a > 



lim z ^l^l ((a, z)) = 0. 



(2.13) 



Proof. We only need to prove (2.13), as the finiteness of |z/|(df) on compact subsets of (0, 00) follows 
from (2.13). It is sufficient to show the claims for 1 > a > 0. We integrate every term on the 
left-hand side of (2.3) from a to 2; and then divide by z. This shows that the limit goes to zero, as 
2; — )■ 00. We first note that 

midx) f°° 
lim / xm(dx) = and lim z~^ / < lim (02;)"^ / m{dx) = 0. 



Hence, 



PZ POO / PZ 

lim / / m(dx)df < lim I / xm(dx) 

J a Jv V J a 



m{dx] 



lim z ^ 



m(da;) 



dt" = and lim z ^ 



X 



m(dx) 



dv = 0. 



So far, we have checked that the terms in (2.3) that do not depend on the tail of the Levy measure 
vanish under the transformation we made, as 2; — )■ oo. Now, we turn our attention to the terms that 
involve the Levy measure of ^. When we'll be dealing with these integrals, the main trick that we 
will use is to change the order of integration. First, we check that 



lim sup z ^ 
< lim sup z' 



lim sup z ^ 



(^In — j m(dx)df 
(-) /. X 



1 



a J V 
z pev 



( In — j m(dx)di; j + limsup 2 (11^ (1) / m{ev,oo)dv 



=(-) 



^(~) / x\ 

lit ( In — ) m(dx)d?7 < limsup 



/ lit ( In — ) df m(dx) 

a Jx/e ^ ^ VJ 



X lim sup z ^ / xm{dx) = 



where we have applied Fubini's Theorem, a change of variables w = \n{x/v) and we have used 
the finiteness of IE[|^i|] and henceforth the finiteness of the quantities Jq Tl^ ^ (^w^ exp{—w)dw and 

tf;'(i). 



=(+) 



Next using Fubini's Theorem and the monotonicity of , we note that for any positive number 



lim sup z ^ 

z—^oo 



Lit I In — ) m(dx)dw 



a Jo 



< lim sup z ^ 

z—^00 

< lim sup z~'^ 

2— >oo 

=(+) 



^0 



X 



:( + ) /, V 



X 



fin — j m(dx)di; = limsup 2; 



z /•ln{z/x) ^ 



={+) 



=(+)/ X 

{w)e'"dw 



xm(dx) + X 



X 

JO 

ln(2/x)Vfe 



Llg [w)e^dw m{d 



X) 



:(+) 



{w)e'"dwm{dx)\ <Il^ {b). 



Since LI (b) decreases to zero as b increases, we see that 



lim z 

2— J-CXD 



a Jo 



^ An — ) m{dx) dv = 0. 



X 



Since rj has a finite mean and m is a finite measure 

rz r={+) 
lim sup z"^ - {v — x)m{dx) dv 

2^00 J a Jo 

_^ /=(+)^ /z\ ri r =(+), , X , \ 

< limsup z n (a) In - I + / - / LI [v — x)m(dx) dv 

( j-z r{x+a)Az 

= lim sup z"^ / m{dx) / LI^ 

2-5-00 \ Jo J aWx 



-(x+a)A2_(+) dy 

[v — x) 

v 



< 



" =(+) 
n, {s)ds 



X lim [az] ^ I m(dx) = 0. 

2— ^-oo 



Similarly, we estimate the following integral 

1 ={+) 
Jo 

/=(+) /z\ poo ^ pw \ 

< limsup 2;"^ I (a) In ( - j + / / — / {w — x)'m{dx) dw dv \ 

pz POO px+a j-^-j 

= limsupz~^ / / / ~n,^ {w — x)dw m{dx) dv 

z^oo J a Jv-a JvVx ^ 

/■"=(+) 1 I 

< / (s)ds xlimsup^;"^ / — m(t> — a, oo) dt> = 0. 

Jo J z^oo Ja V 

As for the remaining two integrals, we split the innermost integrals at the point a; = f + a so that 
{x — v) = (a) and similarly estimate the resulting two terms to get 

limsup^;"^ — (x — f )m(dx) dt> = hmsup ~ l^v i^~'>^)i^i(^^)dwdv = 0. 

Thus, we verify (2.13) and conclude the proof of Lemma 1. □ 

Now that we have established Proposition 1 and Lemma 1, we are ready to complete the proof of 
Theorem 1. 

Proof of Theorem 1. Take an infinitely differentiable function / with compact support in (0, oo) and 
let g{x) := xf'{x). We use (2.6), (2.9), and the identity g{x) = g'{v)dv to get. 



CO 



(j){lD. x)m{dx) = g{x)m{dx) H — - xg'{x)m{dx) 

Jo 2 Jq 

g'{v)Il^ l^ln — j df m(dx) + J J g'{v)Il^ ^In — j dt'm(dx) 



oo \ poo ^ Cp' 



g'{v) \ J m{dx) j dv + J g'{v) ^-^vm{dv) 
+ g\v) (^j^ T[\ ^ {in ^) m(dx)^ dv 



where the interchange of integrals is permitted due to claims of Lemma 1. 
Next, substituting f'{x) = g{x)/x and f"{x) = g'{x)/x — g{x)/x^, we get 



r"£('')/(x)m(dx) = b, r ^m(dx) + ^ T - ^)m{dx 

Jo ^ 2Jo\^^/ 



oo POO / / 



/ I I'q'ix + w) q(x + w)\ =(+) , , 



oo POO / I 



+ / / — - T ^ {w)dwm{dx). 

Jo Jo \ x-w [x-wyj ' 
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Again, using the identity g{x) = g'{v)dv and the fact that 5^ is a function with compact support 
on (0, cxo), we get after careful calculations and an appeal again to Lemma 1 for interchange of 
integration 

Jo Jv ^ ^ \Jo ^ Jo Jv ^ 

poc ^ r=(+) f°° f°° 1 r'=(+) 

+ / g'iv)- / n (v — x)m(dx)dv — / g'(w) / — H (v — x)m(dx)dvdw 
Jo V Jo ' Jo J^ Jo ' 

+ / g'{v)- / {x — v)m{dx)dv — I g\w) I — {x — v)m{dx)dvdw 

Jo ^ Jv Jo Jw ^ Jv 

■.= {9',F,). 

We arrange the above expressions in the form J g'{x)i'{dx), where i^(dx) := Fi{dx) + F2{dx) is the 
same as in Lemma 1. From Lemma 1, we conclude that z/(dx) defines a finite measure on every 
compact subset of (0, 00) and henceforth we consider it as a distribution in Schwartz's sense. Thus 
we get 

POO 

= / £(^)/(x)m(dx) = ig', v) = {g, v') = {xf, v') = (/', xv') = (/, {xv')') , 



for each infinitely differentiable function / with compact support in (0, 00) and derivatives in the sense 
of Schwartz. Therefore using Schwartz theory of distributions for v{dx), we get that xv'{dx) = Cdx 
and therefore 

i^(dx) = (C In X + -D) dx. 

Next, we show that C = D = 0. Note that from (2.13) with a = 1, we have lim^^+oo iy{dv) = 

0. Comparing this with 

0= lim z-^ [ {C\nx + D)dx= lim {Clnz -C + D) 

2— >+oo J^ z—^+00 

we verify that C = D = 0. Thus the proof of Theorem 1 is complete. □ 

The next result is an almost immediate corollary of Theorem 1, and in particular of formula (2.3). 
See also Corollary 3 for a stronger result in a particular case when 77 is a Brownian motion with drift. 

Corollary 1. Assume that condition (2.1) is satisfied. //(t| + > then m(dx) has a continuous 
density on R \ {0}. 

Proof. The absolute continuity of /(^, 77) and boundedness of its derivative on compact subsets of 
(0,00), when cr| + 0"^ > is immediate from (2.3). Let k{x) be the density of m(dx). To show 

the continuity of k{x), we investigate all integral terms in (2.3): all of them, except possibly the 

=(+) =(-) 

ones involving 11^ and 11^ , are clearly continuous. Let us check continuity of these remaining two 
terms. Fix > and v/4 > a > 0. Note that, for any real h such that \h\ < v/4, we have 

j^ln —^j k{x)dx = J j^ln —^j k{x)dx + J j^ln —^j k{x)dx. 

=(+) 

As is continuous and decreasing we verify the dominated convergence theorem applies, as — > 0, 

=(+) 

by bounding in the first term and k{x) in the second. This shows that all integral terms in (2.3) 
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are continuous in v and hence k{v) is continuous. The computation for is the same whereas for 
< we study /(.^, —rf) with the same effect. □ 

3 Exponential functionals with respect to Brownian motion 
with drift 

In the next two sections, we study the special case when rjt = fit + aBf is a Brownian motion with 
drift, so that the exponential functional is now defined as 



oo 



:= / e^'-{pdt + adBt). (3.1) 



We still work under Assumption 1, note that the condition lE[|^7i|] < oo is clearly satisfied. From 
now on, we assume that a > 0, and in order to simplify notations we will write ip{s) = ip^{s). Note 

that formula (3.1) implies J^^o- = o"/^/o-,i, therefore it is sufficient to study the exponential functional 
with a = 1. 

The following three quantities will be very important in what follows 



P 
P 

e 



sup{2 > 
sup{z > 
sup{z > 



E [e^^i] < oo}, 

E [e-"«i] < oo}, (3.2) 
E [e"«^] < 1}. 



In view of (2.2), it is clear that 

p = sup 1^ > : e^^ng(dx) < oo| , p = sup j-z > : e''^n^(-dx) < oo| . 

Thus p > (p > 0) if and only if the measure Il^{dx) has exponentially decaying positive (negative) 
tail. In this case the Levy-Khintchine formula (2.2) implies that the Laplace exponent iplz) can be 
extended analytically in a strip —p < Re(2;) < p. It is clear from (3.2) that < 6 < p. At the same 
time, due to Assumption 1 we have E[^i] = '^'(0) < , which implies that 6* > if and only if p > 0. 
In the next Lemma we collect some simple analytical properties of the Laplace exponent ip{z). 

Lemma 2. Assume that C, satisfies condition (2.1) and that p > 0. Then ip{s) has no zeros in the 
strip < Re(s) < 6. Moreover if ^ has a non-lattice distribution and il){9) = 0, then 9 is the unique 
zero of %Ij{s) in the strip < Re(s) < 6 and the unique real zero in the interval (0,p). 

Proof. Assume that < Re(s) < 9. Since 

we conclude that Re(V'(s)) < V'(Re(s)) < 0, therefore iIj{s) 7^ in the strip < Re(s) < 9. 

Next, assume that il^ifi -\- iy) = Q for some |/ 7^ and ^ has a non-lattice distribution. Then the 
characteristic function of the probability measure e^*'P(.^i G dv) is equal to one at y, therefore it has 
to be a lattice distributed probability measure, see [34, p 306, Theorem 5] which contradicts our 
assumption. 
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In order to prove that 9 is the unique real zero of ip{s) on the interval (0,p), we note that the 
first formula in (2.2) implies that 

^"{s) = a| + I x2e^"ng(dx) > 0, 

therefore '4){s) is convex on (0,p) and it has at most one positive root at 9. □ 
Next, let us introduce two other important objects 

/•oo t2 

:= / e°«*rft, and V := (3.3) 
Jo J2 

We will frequently use the following result, its proof follows immediately from Lemma 2.1 in [27]: 

Proposition 2. Assume that ^ satisfies condition (2.1). For all z & C in the strip —1 < Re(^) < 9/a 
we have E [J^] < 00. 

Our main object of interest is the probability density function of o-? which we will denote by 
k{x) (or by k^^f^{x) if we need to stress dependence on parameters). In the next Lemma, we collect 
some simple properties of k{x). 

Lemma 3. Assume that ^ satisfies condition (2.1). The law of I^^a^ has a continuously differentiable 
density k^^fj{x) which is given by 

k^A^) = / / _=e-^^P(Ji e dy- J2 e dz). (3.4) 
J J (T^/Zirz 



Moreover, both functions k^^^{x) and k'^^{x) are uniformly bounded on R and if fi < then k^^^{x) 
is decreasing on ]R_|_. 

Proof. Expression (3.4) follows by conditioning on ^ and the fact that 

/•oo / /"OO /*oo \ 

J e^W(/xdt + adBt) = N[^j e^^^^ds; J e'^^'^dsj , 
where N{a, b) denotes a normal random variable with mean a and variance b. The continuity of 



T 2 



< 00, see 



k^^tj{x) follows from the Dominated Convergence Theorem and the fact that E 
Proposition 2. 

Next, we observe that the function kHe""" is bounded on M and therefore for some C > we have 



^ ^^\e 5J^P( Ji e dy; J2 G d^) 



< CE[J^''] < 00 



where the last inequality follows from Proposition 2. This shows that we can differentiate the right- 
hand side of (3.4) and obtain 

KA^) = -jj ^^^^'^^^'^^ ^ ^ ^^)' ^3.5) 

and from the above discussion it follows that |A;^^(x)| < CE[J2^^] < 00 for all x G M. Finally, for 
/i < and X > we check that k'^^^{x) < (see (3.5)), therefore k^^„{x) is decreasing. □ 

13 



Our main tool for studying the properties of k^^^{x) will be the Mellin transform of 1^^; which 
is defined for Re(s) = 1 as 

POO 

M^As) ■■= E[(V)^"4{,^.^>o}] = / x'-'k,Ax)dx. (3.6) 

Later we will extend this definition for a wider range of s, but a priori it is not clear why this object 
should be finite for Re(s) 7^ 1. Also, this choice of truncated random variable may seem awkward, 
since we only use the information about the density k^^a{x) for x > 0. However, it is easy to see 
that the Mellin transform A^^,cr(s) uniquely determines k^^„{x) for x > while M._^A^) uniquely 
determines kf^^(j{x) for x < 0. This follows from the simple fact that k^ fj[—x) = k_^^^{x) (clearly 

Iti,a = "1-11,(71 see (3.1)). Moreover, later it will be clear that our definition of the Mellin transform 
is in fact quite natural, since A^^o-(s) satisfies the crucial functional equation (3.13), which will lead 
to a wealth of interesting information about k^^fj{x)- 

As a first step in our study of the Mellin transform M.^A^) obtain its analytic continuation 
into a vertical strip in the complex plane. 

Lemma 4. Assume that ^ satisfies condition (2.1). The function A^^^o-(s) can he extended to an 
analytic function in the strip —1 < Re(s) < 1 + 9, except for a simple pole at s = with residue 
k{0). Moreover, for all s in the strip — 1 < Re(s) < 1 + 6 we have 

M^A^) = ^ + ^ (k{x) - k{0))x'-^dx + k{x)x'-^dx, (3.7) 

and for all s in the strip —1 < Re(s) < it is true that 

1 



00 



Mf,As) = — / x'k'{x)dx. (3.8) 







Proof. First of all, since k{x) is a probability density, it is integrable on [0, 00). Also, due to Lemma 
3, we know that k{x) = k{0) + k'{0)x + o(x) as x — )■ O"*", these two facts imply that Ai^^a{s) exists 
for all s in the strip < Re(s) < 1. 

Next, one can easily check that identity (3.7) is valid for s in the strip < Re(s) < 1. Since 
k{x) — k{0) = k'{0)x + o(x), as X — > 0+ we see that the first integral in the right-hand side of (3.7) 
extends analytically into the larger strip —1 < Re(s) < 1, while the second integral is analytic in 
the half-plane Re(s) < 1. Thus (3.7) provides an analytic continuation of M.^,a{s) into the strip 
— 1 < Re(s) < 1 and it is clear that M.^^o-{s) has a simple pole at s = with residue k{0). 

Next, we note that for —1 < Re(s) < we have 

k{x)x'~^dx= [ {k{x)-k{0))x'-^dx-'^. 



Combining this expression with (3.7) and applying integration by parts we obtain (3.8). 

If ^ = 0, then the proof is finished. However, if ^ > we still have to prove that A^^^o-('S) < 00 
for 1 < s < 1 + 6', and this requires a little bit more work. The proof will be based on certain 
special functions. The confluent hypergeometric function (see section 9.2 in [16] or chapter 6 in [11]) 
is defined as 

n>0 
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where (a)„ = a{a + 1) . . . (a + n — 1) is the Pochhammer symbol. Using the ratio test it is easy to 
see that the series in (3.9) converges for all z G C, thus i-Fi(a, 6, z) is an entire function of z. We will 
also need the parabolic cylinder function, which is defined as 



Dp{z) = 226"^ 



r(¥) 



TT ( p 1 



2' 2' 2 



27r^ 



r(-i) 



1-p 3 ^2 
2 '2' ~2 



(3.10) 



Note that the parabolic cylinder function is analytic function of p and z. See sections 9.24-9.25 in 
[16] for more information on the parabolic cylinder function. We will prove that Ai^^i{s) exists for 
all s in the strip Re(s) G (0, 1 + 6) and everywhere in this strip we have 



Vis) 



-fiVV 



(3.11) 



Let us assume first that Re(s) = 1. Then using (3.4) and (3.6) we conclude that 



E 



X 



V2¥J2 



2 "'2 dx 



-.e 2 



X e ^"'2 dx 



Performing the change of variables x = Uy/J2 and using the following integral identity (formula 
9.241.2 in [16]) 



2 



2 ,2 
uz 



du = r{s)e-D_s{z), Re(s) > 0, 



we obtain equation (3.11). 

Thus, we have established that (3.11) is true for all s on the vertical line Re(s) = 1. Now, we 
will perform analytic continuation into a larger domain. Formulas 9.246 in [16], give us the following 
asymptotic expansions: for z G M 



DJz) 



O (z-'e-- 
0{z-'e-4- 



O ( z'~^ 



e 4 



z — )■ +00, 
Z — !■ — oo. 



(3.12) 



Assume that /i < and s G (0, 1 + 6^) or yU > and s G (0, 1). Then, from (3.12) and the fact 
that Ds{z) is a continuous function of z we find that there exists a constant Ci > such that 

^2 

|e ~^D^s I < Ci for all z > 0. Therefore from (3.11), we conclude that 



\M,,i{s)\ < 



T 2 



and the right-hand side is finite if s G (0, 1 + 6*), see Proposition 2. 

Next, when fi > and s G (1,1 + 0), we again use (3.12) and the fact that Ds{z) is continuous 

in z to conclude that there exists C2 > such that \e~^^D^s (— /^V^) I < C2Z^ for all z > 0. 
Therefore from (3.11), we conclude that 



\M,,,{s)\<C2^E 



'2tx 



2tx 



and the right-hand side is finite if s G (1, 1 + 6*), see Proposition 2. 



□ 
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The next theorem is our first main result in this section. 

Theorem 2. Assume that ^ satisfies condition (2.1) and that 6 > 0. Then for all s such that 
< Re(s) < 9, we have 



M^,^is + 1) + fiMf,,^is) + — (s - l)A^,.,.(s - 1) = 

S Zi 



Proof. Setting 11^ = in (2.3) we find that k{x) satisfies the following integral equation 



(3.13) 



2 Fi{k- v) + h^F2{k- v) + F^{k- v) + Fi{k- v) + ^lF^{k■ v) + —F^{k- v) = 0, 



(3.14) 



where we have defined Fi{k\ v) := k{v) and 

1 

F2{k;v) := - / k{x)dx, F^{k-v) 



1 /■°°=(-) 



n 



Fi{k-v) :-- 
Fe{k-v) :-- 



1 r=(+) 

V 

1 



V 



n. (In 

k{v\ 



V 



X 

k{x 
x"^ 



k{x)dx, F^{k;v) 
dx . 



In j j k{x)dx, 

1 /"^ k{x) 
V ./„ X 



dx. 



(3.15) 



Our plan is to compute the Mellin transform of each term in (3.14). Assume first that 1 < Re(s) < 
1 + min{l, 9}. According to Lemma 4, the Mellin transform of the first term exists in this strip and 
is equal to A^^^o-(s). 

Let us compute the Mellin transform of the second term. We use integration by parts and obtain 
for all 2/ > 



/ v'-^F2{k-v)dv = y'~^ k{x)dx + / 

Jo Jy s-1 Jo 



v" k{v)dv. 



(3.16) 



As ?/ — 7- +00 the first term in the right-hand side of the above equation goes to zero (this follows 
from the fact that the k{x)x'^~^ is absolutely integrable on (0, oo)), thus we conclude that the Mellin 
transform of F2{k] v) is equal to A^^ o-(s) /(s — 1). In exactly the same way one finds that the Mellin 
transform of F^{k] v) is equal to J^^^a{s — l)/(s — 1). 

Let us consider the third term F^{k\ v). Performing the change of variables x ^ yv we find that 



F,{k;v) 



n| \\n{y))kiyv)dy. 



Therefore the Mellin transform of F^lk; v) is given by 



'=(-) 



v'"^F:^{k;v)dv= / (ln(y)) / v'-^k{yv)dvdy 



1 

nj"^ (ln(y))2/-M2/ 



°° =(-1 
m (M)e-(^-i)"dM 



where we have used Fubini's theorem in the first step, performed the change of variables v ^ z/y m 
the second step and applied the change of variables y i— )■ exp(u) in the last step. 
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In exactly the same way we find that the Melhn transform of F^lk; v) is equal to 



Finally, let us consider the sixth term FQ[k\ v). Using integration by parts and the fact that k{x) 
is bounded we find that 



Feik;v) 



k'{x) 



dx. 



X 



(3.17) 



Since k\x) is uniformly bounded on [0, oc) we conclude that F^iji] v) = 0(ln(w)/t>), as f — )■ 0"*", and 
from (3.15) we see that FQ{k]v) = 0{l/v'^), as w — )■ +oo. This shows that the Mellin transform of 
FQ^k] v) exists for 1 < Re(s) < 1 + min{l, 9}. Using (3.17) and integration by parts we find that for 

< f < ^'i < oo 



■"1 yS 

v'-'F,{k;v)dv = ^F,{k;v,) 



Do 



S - 1 



v'~'k'{v)dv. 



(3.18) 



Vo 



From the above discussion we find that the first (second) term in right-hand side of (3.18) converges 
to zero as vi — > +oo {vq — )■ 0"*"), therefore from (3.8) and (3.18) we conclude that for 1 < Re(s) < 
1 + min{l, 6} the Mellin transform of F^lk; v) is given by 



v'-'F6{k;v)dv = - 



s-1 



e _ 2 

v'-'^k'{v)dv = M^,^{s - 2). 



Collecting all the terms in (3.14) we see that for all s in the strip 1 < Re(s) < 1 + min{l, 6} we 
have 



cr, 
'2 



+ 



s - 1 
s-1 



°° =(-) r°° =(+) 

Ul (u)e-(^-i)Mu + / U: (u)e(^-^)"du 



M,,^{s-l) + ^-^M,,^{s-2) 







0. 



Formula (3.13) follows from the above equation by changing variables s H- s — 1 and applying formula 

(2.2) . This ends the proof in the case 6 G (0, 1). If 6' > 1 then (3.13) can be extended from the strip 
< Re(s) < min{l, 6'} to < Re(s) < 6' by analytic continuation. □ 

Remark 1. Note that the functional equation (3.13) is a more general version of the well-known 
functional equation when a = 0, see formula (2.3) in Maulik and Zwart, [27]. Nonetheless, the 
derivation of (3.13) requires the integral equation (2.3) whereas the classical functional equation 

(2.3) in [27] can be obtained by rather simple arguments. 

Theorem 2 will prove crucial for applications. It allows to derive the analytical properties of the 
Mellin transform A^^^o-(s) (such as its behaviour at the singularities and their precise location in the 
complex plane) from the properties of the Laplace exponent ip{s) itself. The next result serves to 
illustrate these ideas. 



Corollary 2. Assume that ^ satisfies condition (2.1) and that 9 > 0. 
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(i) The function A^^o-(s) can be analytically continued into the strip Re(s) G (—1 — p, 1 + p). 
Its only singularities in the strip —1 — p < Re(s) < 1 + 6 are the simple poles at the points 
{-n : < n < 1 + p}. 

(a) If ^ has a non-lattice distribution and 9 < p then M.^^(j{s) has a simple pole at s = 1 + 9 with 
residue 

RiO) ■■= [f^M.AO) + - ^)M,AO - 1)) . (3.19) 

The only other singularities of A4^^„{s) in the strip < Re(s) < 1 + p are poles of the form 
C + n, where n G N and ( is a root of ip{s) in the strip 9 < Re(s) < p. 

(Hi) Consider the "boundary" case 9 = p. Assume that ^ has a non-lattice distribution. If ip{9) < 
then the function A^^^o-(s) is continuous in the strip < Re(s) <l + 9. On the other hand, if 
iIj{9) = and E [i^i exp(6',^i)] < oo, then the function J^fj,,a{s) — R{9)/{s — 1 — 9) is continuous 
in the strip < Re(s) < 1 -\- 9. 

Proof. The proof of parts (i) and (ii) follows easily from Theorem 2 and Lemmas 2 and 4. 

Let us prove (iii). If ip{9) < we use the same argument as in the proof of Lemma 2 and conclude 
that Re(V'(s)) < ?/'(Re(s)) = 'ip{9) < on the line Re(s) = 9; this fact and Lemma 2 imply that 
ijj{s) 7^ in the strip < Re(s) < 9. Since ip{9) < we can use (2.2) and the dominated convergence 
theorem to show that ipls) is continuous in the strip < Re(s) < 9. These two facts and the 
functional equation (3.13) show that A^^,o-(s) is continuous in the strip < Re(s) < 9. 

Finally, let us consider the case when 9 = p and ip{9) = 0. Condition E[^^ exp(6'^i)] < oo and the 
dominated convergence theorem show that the functions ip{s), ip'{s) and ip"{s), which are analytic in 
the strip < Re(s) < 9, can be continuously extended to the right boundary of this strip Re(s) = 9. 
Again, using (2.2) and the dominated convergence theorem one can check that as s — ?• 6' in the strip 
< Re(s) < 9, it is true that 

Note that ip'{9) > due to the convexity of ip{s) on the interval < s < p. Lemma 2 and the fact 
that ^ has non- lattice distribution guarantee that the only zero of ip{s) in the strip < Re(s) < 9 is 
at s = 6'. From here and from (3.20), we see that the function H{s) defined in (3.20) is continuous 
in the strip < Re(s) < 9. 
Let us define 



F{s) = -s (^pM^As) + - l)A^^,.(s - 1)^ 



It is clear from Lemma 4 that F{s) is analytic in some neighbourhood of the line Re(s) = 9, thus the 
function G{s) = {F{s) — F{9))/{s — 9) is also analytic in the neighbourhood of the line Re(s) = 9. 
Next, we use the functional equation (3.13) in the form A^^^o-(s+l) = F{s)/ip{s) and after rearranging 
the terms, we find 

From the above discussion it is clear that the function in the right-hand side is continuous in the 
strip < Re(s) < 9, which ends the proof of part (iii). □ 
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In view of (3.4) it is clear that k^^„{x) depends on the joint distribution of Ji and J2. As we 
will see later in Proposition 3, the Mellin transform A^^,cr(a:) can be expressed in terms of the joint 
moments E[J{' J2]. The next Lemma presents several crucial results on the existence of joint moments 
of this form. Recall that V = Jl/ J2- 

Lemma 5. Assume that ^ satisfies condition (2.1). 



(i) There exists e > such that 



(a) For any {u, s) G in the domain 



E [e'^] < 00. 



(3.21) 



V 



I - 1 < s < 1 + ^, n < o| U |s > 0; n > 0; u < 1 - s| 



we have 



E 



< 00. 



(3.22) 



The function (m, s) G ^ E 



is analytic as long as (Re(s), Re(M)) G T> and it 
is uniformly hounded (Re(s), Re(M)) belongs to a compact subset ofV. 

Proof. Let us prove (i). Denote by Ji{t) = J^e^^ds and J2{t) = /J e^^^ds. It is clear that Ji(0) = 
72(0) = and that both Ji{t) and J2(^) are continuous in t. Since 



dt 







and 



i=0 



dt 



J2{t) 



(3.23) 



t=o 



we conclude that for every x > with probability one, there exists e > such that Ji{t) < xJ2{t) 
for < t < e. This fact and the continuity of Ji{t) and J2{t) imply that 



P (T, <oo)> F{J^ > X J2) 



(3.24) 



where T^ = mi{t > : Jf(t)/J2(t) = x} and as usual we assume that inf{0} = +00. We aim to 
show that for all x > 0, we have g{2x) < g'^{x). 

From the (3.23), we know that Jf(t)/J2(t) — j- as t — )■ 0"*". This fact and the continuity of Ji{t) 
and J2{t) imply that F{Tx = 0) = for all x > and < Ty a.s., for y > x. Using the inequality 
2a2 + 262 > (a + 6)2, we get 

2 (Ji(T,. + t) - Ji(T,.))' + 2J2(T,) > J2(T, + t), 

and we estimate 

^(2x) = P(T2^. < cx)) = P (T,. < cx); 3t > : Ji(T, + t^ = 2xJ2(T, + t)) < 
p(t,. < cx); 3t > : 2 (Ji(T,. + t) - J,{T,)f + 2j2(T,.) = 2xJ2(T, + t))- 

Since JliT^) = xJ2(Tx), we obtain from the above inequality 

g{2x) < P (T, < cx); 3t > : (Ji(T,. + t) - MT,)Y = x{MT, + t)- J2{T^))) = 
F (t^ < oo;3t > : e^^^-J^{t) = xe^^^- J2(t)) = /(x). 
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where J, are the exponential functionals based on = ^r^+t — and we have used the fact that 
the process {C,t}t>o is independent of J-^^. Thus, we have obtained the key inequahty g{2x) < g'^{x). 

Next, let us prove that there exists x* > such that g{x*) < 1. Assume that the converse is 
true, that is g{x) = 1 for all x > 0. In particular, g{n) = 1 for all n > 1. Let An = {3t > : 
Ji{ty = nJ2{t)}. Since = 1 for all n > 1, we conclude that P(n„>iA„) = 1. This implies that 

with probability one there exists a strictly increasing random sequence of positive numbers {tn}n>i 
such that Jf(tn) = nJ2(tn)- Since J2(^n) is an increasing sequence, we conclude that as n — >■ +oo 
we have P(Ji(t„)^ — )■ +oo) = 1, and due to the fact that Ji(t„) < Ji, we arrive at a contradiction 
P(Ji = oo) = 1. 

Thus, we have proved that there exists x* > such that g{x*) < 1. For x > x*, let us define 
> to be the unique integer number such that 2^ < x/x* < 2^+^. Applying the inequality 



g{x) < g{x/2) exactly A^ times, we obtain g{x) < g{x/2 



Using the fact that g{x) is a decreasing 



function and that x* < x/2 , we conclude that g{x) < g{x*) , and since x/{2x*) < 2 , we see that 
for all X > X* 

g{x) < e-'=^ 

where c = -\\i{g{x*)) /{2x*) > 0. This fact and (3.24) imply that F{V > x) < exp(-cx) for all 
X > X*, thus (3.21) is true for any e G (0, c). This ends the proof of part (i). 

Let us prove (ii). Assume first that u < 0. Then using Holder inequality we get 



E 



i(u+s-l)' 



E 



< (E [V-^P]y (E 



T2V 
'^2 



From part (i), we know that V has finite positive moments of all orders. Then it suffices to choose 
g = 2(1 - s)~^ for -1 < s < 0, g = 2 for < s < 1 and g = | + 



2(s-l) 



for 1 < s < 1 + 6' to conclude 



that (3.22) holds. 

Assume next 0<s<l, m>0 and u < 1 — s. Then with p = and g = (1 — u)~^ we have 



E 



T-« A{U + S-1) 



J. 



-(-IH — ^1 



1-u 



< 00, 



□ 



due to Proposition 2 and the fact that (—1 + G (—1, 0]. 

Now we are ready to present several integral expressions for the Mellin transform A^^^o-(s). These 
expressions are interesting in their own right, but they will also lead to an important result about 
the exponential decay of A^^,o-(s) as Im(s) — > 00 (Theorem 3 below). Note that due to the identity 

o"-^^j/cr,i we have A/l^^o-(s) = cr*^^A/l^/o-,i(s), therefore it is enough to state the results for cr = 1. 



Proposition 3. Assume that C, satisfies condition (2.1). 



(i) For -1 < Re(s) < 1 + ^ 



2-l(^+i)r(g) 
r (1(^ + 1)) 



E 



(ii) For /i < and -1 < Re(s) < 1 + 9 

2-1(^+1) 



27ri 



^2 



(s-1) 



(3.25) 



^iMF(i(n + s + l)) 



•E 



(2/i2)-tdM. (3.26) 
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(iii) For /i > and -1 < Re(s) < 1 + ^ 

-M.^,,{s). (3.27) 

The proof of this Proposition is quite technical, therefore we have divided it into several steps. 
First of all, in Lemmas 6, 7, 8, we establish several technical results which will be needed in the proof 
of Proposition 3, and also useful later. 

Lemma 6. 



2 2 2 



(i) For every e > and a < b there exists C = C{e,a,b) > such that 

\T{x + iy)\ < Ce-(5-^)l^l 

for all a < X < b and \y\ > 1. 
(a) For every e > and a <b there exists C = C{e,a,b) > such that 

\r{x + iy) \ > Ce~(i+')l^l 

for all a < X < b and ?/ G M. 
Proof. We start with the following asymptotic expression 

1 



r(a; + iy) 



27r\y\ 



y — 7- OO 



(3.28) 



which holds uniformly in x on compact subsets of M, see formula 8.328.1 in [16]. Part (i) follows 
easily from (3.28) and for part (ii), we use the additional fact that r(s) has no zeros in the entire 
complex plane. □ 

Lemma 7. For fi<0, Re{w) < \ and < Re(s) < 1 - 2Re{w) 



(x-iJ.y) 



2r^x'-^z'"-^dxdz = 2"'~^(-/iy) 



\2w+s 



^ V27ri 1(1-^) 
Proof. We change the variable of integration z ^ and find that for a > and Re(tf) < 1/2 



(3.29) 



^.au^-\-w^^ = a'^-'^vi^ -w\. 



Then for 2Re(ii;) + Re(s) — 1 < we can apply the Fubini's theorem and obtain 

(x - /iy)'"'^ V-Mx 



X {x — fiy)'^ -l -| 

; 2i x''' z"^' dxdz 



TTZ 



^2^-'^t(- -w 
/2^ V2 



2wA 



-1 ^ =25-™r 
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a/2^ 



w 



w 



{x + lf'"-^x'~^dx 

r(s)r(l-2M;-s) 
r(l -2m;) ' 



where in the last step we have used the beta-integral identity (see equation 3.194.3 [16]). Formula 
(3.29) can be derived from the above equation by application of the Legendre duplication formula 
for the gamma function (see formula 8.335.1 in [16]). □ 
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Lemma 8. Assume that ag < ai and 6 G C are such that Re(6) G (0,1) U (1,2). Recall that 
iFi{a, b, z) denotes the confluent hypergeometric function defined by (3.9). For each e > there exist 
a constant C = C{ao,ai,b,e) > and a constant D = D{ao,ai,e) G (0, |) such that for all a E C 
with Oq < Re(a) < Oi and all z > 



\iFi{a,b,-z) \ < Ce^-+^'^'"('^)L 
Proof. We start with the following integral representation 



iFi(a,6, -z) = — r(6)z^-V 



e-^-^d _ yj- 



"dw, 



(3.30) 



(3.31) 



7+iR 



which holds for z > 0, Re(6) > and 7 > 1. This representation follows from formula (7) on page 
273 in [11] and the identity iFi{a,b, —z) = exp{—z)iFi{b — a,b,z) (see formula (7) on page 253 in 

[11])- 

Next fix e > and assume that Re(6) G (1,2) and z > 1. We also denote 7 = 1 + e. Then 
changing variables w h- )■ 7 + it we obtain from (3.31) 



iFi{a,b,z)\ 



1 



:T{b)z^-''e'' 



27n 
< Ci{b)e 



°° / \ a-b 

itz/^, I ■±\~b 



e**'(7 + zt)- 



(7 + zt)-' 



dt 



(7 + zt)-^ X (l-(7 + it)-i) 



dt. 



(3.32) 



Note that the set {(7 + it)"^ : t G M} C C is a circle with centre (27) ^ and radius (27) ^. 
Therefore the set {1 — (7 + it)^^ : t G M} C C is a circle with centre 1 — (27)""'^ and radius (27)""^. 
Recall that 7 = 1 + e > 0, therefore this last circle does not touch the vertical line iM and we have 



D = max I |arg(l - (7 + zt) ^)| | < 
At the same time, we have for alH G M 



(3.33) 



1 + e 



< 1 1 - (7 + it) 



-1 1 



g^ + t^ 
72 + ^2 



< 1. 



The above two estimates and the equality \u'"\ = |ii|i^'=(*')el^''s(")^^™('')l (which is valid for all m G C and 
f G C with Re(M) > 0, Re(w) > 0) show that for all t G M, we have 



< 6*2(00, ai, 6, e)e 



1 + e 



ao— Re(b) 



{i-{^+ttr') 

where 

^2(00, 0,1, b, e) = max 1 
Using (3.32) and (3.34), we conclude that 

\,F^ia,b,-z)\ < Ci(6)C2(ao,ai,6,e)e^^+^l^™('^ 
= C(ao,ai,6,e)e"+^l^"^('^)l. 



D\Im{a)\ 



(3.34) 



1 + e 



ai — Rc(6) 



|(7 + 2t)"^|dt 
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Note that the integral appearing in the above estimate converges since Re(6) G (1,2). This proves 
(3.30) for z>l. 

Assume next that z G (0, 1). Using (3.31) with 7 = {1 + e)/ z and changing variables in the 
integral w ^ (1 + e + it)/z we get 

,F,(a, i, = q^r(„^ (1 + . + it)-' (1 - iTTTa)""*- 

Now we can proceed as in the case when z > 1 noting that the set {1 — z{l + e + it)~^ : t G M} C C 
is a circle with centre 1 — z{2{l + e))~^ and radius ^(2(1 + e))~^. As < z < 1 one can see that this 
only improves all the estimates above. For example, the estimate (3.33) also holds true and for all 
t G M, we have 



e l + e-z . ,^ . , 1, /(l + e-z)2 + t2 

< — < l-2(l + e + it)-H = J^- . < 1. 



1 + e 1 + e -I ' 'I V (l + e)2 + t2 

Therefore (3.30) is also true for z G (0, 1). 

Finally, we consider the case when Re(fe) G (0, 1). One can see that this case follows easily from 
the already established result valid for Re(6) G (1,2) and the following identity for the confluent 
hypergeometric function 

h\F\(a, b, —z) = aiFi{a + l,b + 1, —z) + {h — a)iFi(a, 6 + 1, —z), 

see formula 9.212.3 in [16]. □ 
Proof of Proposition 3. The equation (3.25) follows from (3.11) and the fact that 

v^2-f 



I?-.(0) 



r(i(. + i)) 



see the definition of the parabolic cylinder function (3.10). 

Let us prove (ii). Assume first that s G C is a fixed number which satisfies Re(s) G |). Using 
Lemma 3 and Fubini Theorem we find that the Mellin transform of k{x) is given by 

M^,i{s) = jj F{s, y, z)P( Ji G dy, h G dz), (3.35) 



where 

F{s,y,z) = —= I x^~^e~ ' 2^ dx. 

According to Lemma 7, since Re(s) G (^, |) the Mellin transform of F(s, 2) in the ^-variable exists 
for all w such that < Re(w) < 1/8 and is given by 

G{s, y, w) := r F(., y, z)z-'dz = 2-^(-/iy)^-+-^ ^^^^^^/ ~ ~ (3.36) 
Jo ^ ~ "^J 
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Using Lemma 6, we find that for every s sucli tliat 1/4 < Re(s) < 3/4 tliere exists C = C{s) > 
sucli tliat for all w with Re(u;) = 1/16 we have 



\G{s,y,w)\ < C|y|^'=W-ie-l^™('")l. 



(3.37) 



Therefore as \G{s,y,w)\ is absolutely integrable along the line w = ^ + iM then F{s,y,z) can be 
written as an inverse Mellin transform 



F{s,y,z) = ^[ G{s,y,w)z ""dw. 
2-111 y 1 +iM 



From the above identity and (3.35) we find that 



1 
27ri 



OO POO 



y=0 Jz=0 JtoG-y+i: 



G{s,y,w)z~'"dwF{Ji e dy,J2 G dz). 



Due to (3.37) and the fact that E,[jf'^^^^~'^^^ J2^^^^] < oo since 1/4 < Re(s) < 3/4 (see Lemma 5), 
we conclude that the function G{s,y,w)z~^ is absolutely integrable with respect to the measure 
dw X P( Ji G dy, J2 G dz). Thus, we can apply Fubini's Theorem to the right-hand side of the above 
equation and with the help of (3.36), we obtain 



27ri 



r(l-2u;-s) 



rri 



w) 



Next, we perform a change of variables w ^ ^{1 
obtain from (3.38) 



E [j^'^+^-Vs""] 2"'-i(-/i)'"'+^-Mw. (3.38) 
s) (recall that s is a fixed number) and 



u 



2-5(«+i)r(s) 

2^ri 



Ffn) 



,+iMF(i(M + S + l)) 



E 



J-U J 



;{U+S-1) 



{2fi^)-^du. 



(3.39) 



(u+s-l) 



IS 



For s fixed, such that 1/4 < Re(s) < 3/4, we know from (ii) Lemma 5 that E "J, 

a bounded analytic function of u everywhere in the strip —1 < Re(M) < 17/18 — Re(s) and hence 
bounded on Re(M) = 7/8 — Re(s). The ratio of Gamma functions F(u)/F(|(m + s + 1)) decays 
exponentially (and uniformly) as Im(u) — )■ cxd in the strip —1 < Re(u) < 17/18 — Re(s), and it has 
a unique simple pole at m = 0, coming from T{u). Thus we can shift the contour of integration in 
(3.39) 7/8 — s + iM I— )■ — 1/2 + iM and taking into account the residue at m = we finally obtain 



+ 



2-|(^+i)F(s) 

r (1(^ + 1)) 

2-|(s+i) 



E 



Or. 



(3.40) 



T{s)T{u) 



27ri J.^,^T{\{u + s + l)) 



E 



J-U T2("+«-l) 



(2/i2)-tdu. 



According to Lemma 5, E 



J-U jZ 



{u+s-l) 



is a bounded analytic function for Re(w) 



-1/2 and 



— 1 + e < Re(s) < 1 + ^ — e, for any e > 0. Due to Lemma 6, the ratio of Gamma functions 
F(u)/F(|(n + s + 1)) decays exponentially as Im(n) — 00, Re(u) = —1/2 and uniformly in s if 

— 1 + e < Re(s) < 1 + 9 — e. Therefore, the right-hand side in (3.40) defines a meromorphic function 
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in the strip —1 < Re(s) < 1 + 6, which has a unique simple pole at s = (which comes from 
r(s)), and we can apply analytic continuation and conclude that (3.40) is valid for all s in the strip 
— 1 < Re(s) < 1 + 9. This ends the proof of part (ii). 

Finally, let us prove (iii). Assume first that < Re(s) < 1. We use formulae (3.10) and (3.11) to 
find that 



M.i(s) + M^ 



Vis] 



2-i(s-i) 



E 



(3.41) 



From the above formula and the identity e~^iFi(a, b, z) = iFi(6 — a, b, —z) (see formula (7) on page 
253 in [11]), we conclude that (3.27) holds true for < Re(s) < 1. Now our goal is to check that 
formula (3.27) can be extended into the wider strip —1 < Re(s) < \ + 6. 

Assume that 5 > is a small number and that —l+5< Re(s) < 1 + 9 — 5. It is clear that 
we can find p = p{6) > 1 such that for all s in the strip —l + 6< Re(s) < 1 + 6 — 6, we have 
(Re(s) — l)pE (—2,6'). Define q = p/ip — l)- According to Lemma 5, we can find e > small enough 



such that E exp{eq^V) 
C = C{6)>0 such that 



< oo. Using Lemma 8, we see that there exists D = D{6) G (0,7r/2) and 
br all s in the strip —l + d< Re(s) < 1 + 6' — 5, we have 



li^i 



1 - s 1 



1^' 



'2' 



V 



Therefore, we can use Holder inequality with p and q defined as above and estimate the expectation 
in the right-hand side of (3.27) as follows 



E 



Jo l-Tl 



1-S 1 



2 2 2 



< 



^gD|Im( 



E 



■ l(Re(s)-l)p- 
•^2 



E 



< OO, (3.42) 



where in the last step we have used the fact that |(Re(s) — l)p E (—1,9/2). This shows that the 
expectation in the right-hand side of (3.27) is well-defined for all s such that —1+6 < Re(s) < 1+9—6, 
and since 5 > is an arbitrary small number, we can extend the validity of this equation into the 
whole strip -1 < Re(s) < 1 + ^. □ 



The next theorem is our second main result in this section and it opens the way for the application 
of powerful complex-analytical tools. 

Theorem 3. Assume that ^ satisfies condition (2.1). For any /i G M and any small number 6 > 0, 
there exist constants A = A{fi, a,6) > and B = B{fi, a,6) > such that 



\M,As)\<Ae-''\'"'^'^\, 
for all s eC such that Re(s) G (-1 + 6,1 + 9-6) and |Im(s)| > 1. 



(3.43) 



Proof. Note that ha = o-Ie i, hence M-aais) = cr'^ M.e As), therefore without loss of generality 
we can assume a = 1. 

Since —1/2 is not a pole for r(s), we use Lemma 6 and conclude that there exists Ci > such 
that for all ?/ G M 



r(- 



\y\ 



(3.44) 
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At the same time, from Lemma 6 we find that for all x G (—1,1 + ^) and y G M there exists C2 > 
such that 



r(x + iy) 



(3.45) 



E 



First let us assume that /i = 0. Then (3.43) follows immediately from Lemma 6 and (3.25) since 
J2 ^ < C((5), for Re(s) G { — 1 + 5,1 + 6 — 5). The latter is obvious from Proposition 2 for J2. 

Next, assume that < 0. Thanks to (3.26) and Lemma 5, we get that 

T{s)\ f \V{-\ + iy) 



\M,,M\<C{^^,5,e) 



\n\{s + l))\ 



Wis) 



|r(HI + i2/ + ^)) 



-dy 



(3.46) 



From Lemma 6, we deduce that for —1 < Re(s) < 1 + 6* and |Im(s)| > 1 there exists C3 > such 
that 



Wis) 



W{\{s + 1))\ 



which shows that the first term in (3.46) is decaying exponentially as Im(s) -+ 00. Next, from Lemma 
6, we know that for —1 < Re(s) <1 + 9 and |Im(s)| > 1 there exists C4 > such that 

|r(s)| <C4e-i^i^"^(^)i. 

Using this fact and estimates (3.44) and (3.45), we see that for |Im(s)| > 1 



Wis) 



1^(2 (2 +12/ + *)) I ^2 



e -^lyl+ily+imWI 



dy 
dy 



e 32 \y\dy = C^e~ 



|Im(s)| 



The above estimate shows that the second term in (3.46) is decaying exponentially as Im(s) -+ 00, 
which ends the proof in the case /i < 0. 

Finally, let us consider the case when /i > 0. In the proof of part (iii) of Proposition 3 (see 
inequality (3.42)), we have established that for every 5 > there exist constants D = D{5) G (0, vr/2) 
and C = C{5) > such that for all s in the strip —l + 5< Re(s) < 1 + 6' — 5 we have 



E 



J. 



1 _^ 

2 '2' 2 



< 



Ce^|Im(f)|_ 



Then from (3.27), we find that for all s in the strip —1 + 5 < Re(s) < 1 + 6 — 5 

Vis) 



K^M,i(s) < 



n\{s + i)) 



(3.47) 
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Using Lemma 6 and the fact that D < 7r/2, we conclude that there exist Ci > such that for all s 
in the strip —l + 6< Re(s) < 1 + 6 — 6 and |Im(s)| > 1, we have 



r(i(^ + i)) 



Therefore the first term in (3.47) can be bounded by 



and it decays exponentially as Im(s) — oo since 7t/2 — D > 0. This ends the proof in the case > 0, 
since we have already established that the second term in (3.47) decays exponentially to zero. □ 

Corollary 3. Assume that satisfies condition (2.1). The function k{x) is infinitely differentiable 
on R \ {0} and k{x) E C^(M). 

Proof. The fact that k{x) G C"'^(]R) was already established in Lemma 3. Assume that x E (0, oo). 
Applying Mellin transform inversion, we find that 



where the integral converges absolutely since (3.43) guarantees exponential decay of Ain^ais) on the 
line 1/2 + iR. This exponential decay also guarantees that for every n > the functions 



are absolutely integrable along the line 1/2 + iIR, which shows by differentiation under the integral in 
(3.48) that k^^aix) E C°°(0, oo). Noting that -J^,^ = /_/,,^ we deduce that k^^^ix) E C"^(-oo, 0). □ 

Corollary 4. Assume that ^ satisfies condition (2.1) and that p > 0, 6 > 0. For any /i G M and any 
small number 6 > the estimate (3.43) holds uniformly in the strip Re(s) G ( — 1 — p + 6,1 + 6 — S). 

Proof. The statement about the exponential decay follows from Theorem 3, the functional equation 
(3.13) and the fact that ipi^) = O^z"^) uniformly in the strip Re(s) G (—1 — p + 6,1 + 9 — 6). The 
latter fact follows from (2.2) (see also Proposition 2 in [3]). □ 

Theorem 3 is very important for several reasons. First of all, as we have seen in Corollary 3, it 
implies smoothness of k{x) on R \ {0}. This should be compared with the case a = 0, where it is 
known that k{x) may be non-smooth on (0, oo). For example, if ^ has bounded variation and negative 
linear drift p^, then k{x) may be non-smooth at point —1/p^, see Proposition 2.1 in [9] and remark 
2 in [19]. Secondly, as we will see in the next section. Theorem 3 together with Theorem 2 will allow 
us to use simple techniques from Complex Analysis, such as shifting the contour of integration in 
the inverse Mellin transform, to prove rather strong results about the asymptotic behaviour of k{x) 
as x — )■ O"*" or X — )■ +oo. 




(3.48) 



n-l 



Ylis + i)M^,,ais) 



i=0 
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4 Applications 



In this section, we present several applications of the results obtained in the previous section. We are 
still working under the same assumptions as in Section 3, i.e. we consider the exponential functional 
Ifj^^a defined by (3.1) under the assumptions: IE[|^i|] < oo, E[,^i] < and cr = a,, > 0. 

Our main tools are the meromorphic extension of A^^,^, Tauberian theorems and Mellin inversion 
with shifting of the contour of integration. We will also use the functional equation (3.13) and the 
estimate (3.43) developed in Section 3. In Theorem 4 we derive some asymptotic results for k{x) 
as X — 0, while in Theorem 5 we discuss the behaviour of IP(/^,o- > x) and k{x) as x — )■ oo, thus 
strengthening significantly some of the results of Lindner and Mailer [25, Theorem 4.5] in this special 
case when rjs = fis + aBg. We note that under further assumptions much stronger results are within 
reach for the asymptotic behaviour of IP(/^,o- > x) and k{x) both as x — ?■ and x — ?■ oo. In order 
to illustrate the techniques, we choose to work with a rather simple but nevertheless very useful for 
applications class of processes ^ which have hyper-exponential jumps (see [7, 8, 20]). The same results 
can be easily generalized to more general class of Levy processes with jumps of rational transform 
(see [19]). 

Finally we point out that IP(/^,o- > x) can be associated to ruin probability for certain actuarial 
models, see for example Theorem 4 in [1]. 

4.1 General results about asymptotic behaviour of k{x) 

Our first theorem in this section deals with the asymptotic behaviour of k{x) at zero. As usual, we 
define the "floor" (or "integer part") function as \_x\ = max{n G Z : n < x}. We recall that p is 
defined by (3.2): if p > then the Levy measure of has exponentially decaying negative tail with 
the rate of decay equal to p. 

Theorem 4. Assume that ^ satisfies condition (2.1) and that 6 > 0, p > 0. Then for every integer 
m > and e G (0, 1) such that m + e<l+pwe have 

m , 

(x) = V4x" + 0(x™+n, asx^O, (4.1) 

where the coefficients {bn}n>o ore defined recursively: 6„i = 0, bo = k{0) and 

2 

bn+i = ifJ'bn - 'ip{-n)bn-i) , for < n < p. (4.2) 

In particular, k{x) G C^'^^^^(K), and if p = oo then k{x) G C°°(M). Moreover, as Remark 5 shows 
k"{0) may fail to exist. 

Proof. Recall from Corollary 3 that Ai^^f^{s) is analytic in the strip Re(s) G (—1 — p,l + 6) and has 
simple poles at all negative integers —n such that < n < 1 + p. Define 

(^n = Cin{p, cr) = Res(A^^^o-(s) : s = — n), < n < 1 + p. (4.3) 

Choose c < 1 + p, such that c ^ N. We start from the Mellin transform inversion formula (3.48), use 
the fact that A^,i,o-(s) decays exponentially as Im(s) — > oo (and uniformly in Re(s)) and shift the 
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contour of integration 1/2 + iR i— >■ — c + iM while taking into account the residues at points —n to 
obtain 



27ri y^=i , iK ivri y_c+iR 



(4.4) 



0<n<c 



The integral term in the right-hand side of the above equation can be estimated as follows 



J~c+m Jr 



dt, 



therefore this term is 0{x'^) as x — ?■ 0+. 

Let us derive a recurrence relation for the coefficients a„. First of all, from Lemma 4 we find that 
ao = k{0) (this fact is also obvious from (4.4)). Next, from the definition (4.3) and the fact that all 
the poles are simple, we find that 



s + n 



0(1), s^-n. 



(4.5) 



Using formula (4.5) and the functional equation (3.13), we find that as s — we have 



A^M,<x(s + l) + /"- + 0(l) + 



a 



ai 



0(1) = 0. 



Due to the fact that ■ip{s)/s — )■ ^''(O) = < oo, as s — j- 0, and that A^^_o-(s-|- 1) — )■ 1, as s — )■ 0, we 
conclude that fiao — = 0. Following the same steps and considering the functional equation 

(3.13), as s — 7- —n, we find that the coefficients satisfy the recurrence relation 



tfj^—n) 



_^ -a„,„i + fian + ~ = 0' > 1- 

Therefore, if we define 6„ = bn{fJ',cr) = n\an{fi,a) then from the above equation, we obtain the 
recurrence relation (4.2). 

Combining all the above results we see that we have established (4.1), but only in the one-sided 

sense x — )■ 0^. Using the fact that I^^a = ~I-fj.,a and repeating the above arguments, we obtain 
k^A-x) = k-^A^) = + 0(x™+^), as 0+. 



n=0 



Clearly, (4.1) would be true if bn{—fJ',o') = (— l)"6„(/i, a). This fact can be easily verified: using the 
recurrence relation (4.2), we check that 6,1 (/i, a) is a polynomial in of degree n, which is odd (even) 
if n is an odd (even) number. This ends the proof of asymptotic formula (4.1). 

Finally, formula (4.1) and the fact that A;(x) G C^{R \ {0}) imply that k{x) E 0^+LpJ(m), which 
ends the proof of Theorem 4. □ 

Remark 2. To the best of our knowledge this is the first general result on the behaviour of k^ Ax) as 
X — in the case cr > 0. At the same time there are several recent results concerning such behaviour 
when a = 0, see [19, 21, 28]. 
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Note that if the process ^ is spectrally positive, or more generally, if n^(dx) restricted to (—00, 0) 
has exponential moments of arbitrary order, then k{x) G C°°(]R). 

Our next result provides an extensive account of the asymptotic behaviour of f{I^,a > x) as 

X —7- +00. 

Theorem 5. Assume that ^ satisfies condition (2.1) and that 6 > and ^ has a non-lattice distri- 
bution. 

(i) If one of the following conditions is satisfied 
(a) 9<p, 

(h) e = p, ipie) = and E[^lexp{e^i)] < 00, 
then 

P(/m,<x > x) = Cx-^ + [x-^] , X +00, (4.6) 
where C = —R{9)/9 and R{9) is defined in (3.19). 
(a) If9 = p and iIj{9) < then 

P(/^,^ >x) = o (x"^) , X +00. (4.7) 

Moreover, if p < then the asymptotic expressions (4.6) and (4.7) can be differentiated and leads to 
an asymptotic expression for k{x). 

Proof. Let us prove part (i). For x > we define 

X 

F{x) := j P(/^,, > y^)dy. 


Using integration by parts in the same way as we did above when dealing with the Mellin transform 
of the F2{k; v) term in the proof of Theorem 2 (see also equation (3.16)), we find that for all s in the 
strip Re(s) G (0, |) 



POO 

F{s) := / x'-^dF{x) = s~^7W^,,(l + 29s) 
Jo 



(4.8) 



Due to Corollary 2, the function F{s) — C/{l/2 — s) is continuous in the strip < Re(s) < 1/2, 
therefore we can apply Wiener-Ikehara Theorem (see Theorem 7.3 in [2]) and conclude that as 

X —7- 

F{x) = 2Cy/x + O (y/x) . 

Using the above asymptotic expression, the fact that P(/^,o- > y^) is a decreasing function of y and 
applying the Monotone Density Theorem we obtain (4.6). 

The proof of part (ii) is very similar: now we use Corollary 2 to find that F{s) is continuous 
in the strip < Re(s) < 1/2, therefore by applying Wiener-Ikehara Theorem we conclude that as 
X —7- +00, we have F{x) = o{^/x), and applying the Monotone Density Theorem we obtain (4.7). 

If p < then from Lemma 3 we know that k{x) is a decreasing function on (0, 00), therefore we 
can apply the Monotone Density Theorem to (4.6) or (4.7) and obtain the corresponding asymptotic 
expression for k[x). □ 
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Remark 3. Note that, for allt > 0, 

Jo 

where /^^ = /^^o- and /^^ independent of ^e^% e^'dr/^y Recalhng that in [25, Prop. 4.1 ; Theorem 

4.5] the authors use — ^ for ^ in the definition of J^^o-, we point out that the authors supplement the 
theory of random recurrent equations developed in [15] to deduce general results for the behaviour 
of P (/(^, f]) > x). For the case when r]t = + aBt their result translates to 

lim x^P >x) = C+>0; lim < -x) = > 0; hm x^P > x) = C+ + C_ > 

x— >-oo ' x—^oo ' x—¥oo ' 

under the conditions (following our notation) that either p>0>lor^<l<p. Our assumptions 
are much weaker, see (i) Theorem 5 and we also compute the constants C±. 

Moreover, assuming /i < 0, or otherwise working with = —Ifj,.a, one can show that —R{0) > 
0, thus C+ > and hence C+ + C_ > 0. To prove that —R{6) > 0, we consider two cases: when 
^ > 1 this follows directly from (3.19) and the fact that ip'{6) > 0, and when 6 G (0, 1) this follows 
from (3.19) and the fact that {6 — l)M^,aiO — 1) > (the latter is true due to (3.8) and the fact 
that k{x) is decreasing, see Lemma 3). 

Finally we note that, despite dealing with the asymptotics of t]) for general t], the methodology 
in [25] cannot seemingly be improved to yield stronger results for the special case when 77 is a Brownian 
motion with drift. 

Remark 4. The case when a = has been completely dealt with in [27, 31, 32]. We note 
that the technique applied there again relies on the random recurrence equations studied in [15] 
and the authors are able to obtain results in part (i), condition (b) under the weaker assumption 
E [^1 exp(6'^i)] < 00. A recent paper by V. Rivero [33] addresses the case when the process ^ has 
convolution equivalent Levy measure, the main tools are fluctuation theory of Levy processes and an 
explicit path-wise representation of the exponential functional. 

4.2 Case study: processes with hyper-exponential jumps 

In this section, we will show how our methods can be extended to derive quite strong results about 
the density of the exponential functional, provided that we impose additional restrictions on the 
Levy process ^. In particular, we will need more information about the analytical structure of the 
Laplace exponent iIj{z). Our purpose in this section is not to prove the most general results possible, 
but rather to present the ideas and give the flavour of the results which can be derived. 

Let us consider a simple (but very useful) class of processes having hyper-exponential jumps (see 
[7, 8, 20]). In this case the Levy measure of a process ^ is essentially a mixture of exponential 
distributions 

N N 

n^(dx) = l{^>o} J2 «ne-^"^dx + l|^<o} ^ a^e'^-Mx. (4.9) 

n=l n=l 

where all the constants a„, dn, Pn, f>n are strictly positive. Since A = ng(M) < 0, the process ^ can be 
represented as Brownian motion with drift plus a compound Poisson process 

N{\t) 
n=l 
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where N{t) is the standard Poisson process and Yi are i.i.d. random variables with distribution 
P(li G dy) = X~^Il^{dy) . Note that fi^ is the hnear drift of the process ^, it is easy to relate it to the 
constant = E[^i] by 6^ = /i^ + AE[Yi]. 

Formula (4.9) implies that the Laplace exponent of a hyper-exponential process is a rational 
function of the form 

2 N N ^ 

(4.10) 

n=l ^ ' n=l ^ ' 

For hyper-exponential processes, it is known (see [7, 20]) that the equation ip{z) = has only 
real simple solutions. Denote the positive solutions as {Cm}i<m<M, where we assume that they are 
arranged in increasing order. It is also known that M = iV -f 1 if (i) cr^ > or (ii) = and /i^ > 0, 
and M = N otherwise (see [19]). Note that in our previous notation (3.2), we have = (i, p = pi 
and p = pi. 

Using this information about zeros and poles of ip{z) and the functional equation (3.13) it is easy 
to see that Ai^,a{s) can be extended to a meromorphic function, with poles at the points 

{Cm + n : m > 1, 1 < n < N} U {-p„ - m : m>l, l<n<iV}U {-m : m > 0}. 

If we further assume that 

Q - ^ Z for all 1 < 2, J < M and z ^ j, 
Pi ^ N and pi — pj ^ Z for all 1 < i, j < N and i ^ j, 

then it is clear that all the poles of are simple. Let us introduce the following notations 

Cij = --^Res (7W^,a(s) : s = j + (i), l<i<M, j> 1, 

bij = {1 + pi)jRes{M^,a{s) : s = -J - Pi) , l<i<N, j > 1, 

recall that (a)„ = a(a+l) . . . (a+n — l) denotes the Pochhammer symbol. Our next goal is to compute 
coefficients Cij and bij in terms of the Mellin transform A^^ o-(s). Let us fix i such that 1 < i < M. 
Since Q is a simple root of a rational function ip{z), we have ip{z) = ■ip'{Q){z — Q) + 0{{z — QY) as 
z — > Q. This fact and the functional equation (3.13) show that 

c,i = ^ [fiM,AQ) + y(0 - i)M,AQ - 1)) • (4.12) 

Next, using the functional equation (3.13) and the same technique as in the proof of Theorem 4, we 
obtain a recursion equation 



1 f a 



^(j + 0) V" 2 



2 



/^Cij + ^Qj-i , j > 1, (4.13) 



where we have defined q^o = 0. Next, let us fix i such that 1 < i < N. Formula (4.10) implies 
that ip{z) has a simple pole ai z = —pi with residue Oj. Again, we use this fact and the functional 
equation (3.13) to conclude that 

ki = -^?^-^M,<x(l-A)' (4-14) 

pi 
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and 



2 

= ^ ^f^^i'j ~ ^(^-^ " Pi)Kj-l) ' J > 1, (4-15) 



where we have defined 6j q = 0. Recall that the coefficients bj = j!Res(A^^^o-('5) : s = —j) can be 
computed via the recurrence relation (4.2) 

Our main result in this section is the following Theorem, which provides a complete asymptotic 
expansion for k{x) as x — ?■ 0^ and at x — ?■ +oo. The corresponding expansions as a; — )■ 0~ and 
X — — oo can be obtained by considering k^^^^^x) = k^^^^—x). 

Theorem 6. Assume that C, is a process with hyper- exponential jumps (4.9) , which satisfies < 
and for which conditions (4.11) are satisfied. Then for every c > 

bi . , x^+P^ 



0<j<c ■'' i=l j>l ^ P"^'^ 

^^'^{j+C.<c}Cij^^^ + O [x-') , as x^+oo. 



(4.16) 



Proof. From (4.10) it is clear that ip{z) = 0{z^) and l/ip^z) = 0(1) as Im(2) — )■ oo, |Im(2)| > 1, and 
that these estimates are uniform in Re(s). Therefore, using Theorem 3 and the functional equation 
(3.13) we see that A^^,o-(s) decays exponentially as Im(s) — j- oo, Im(s) > 1, and uniformly if Re(s) 
belongs to a compact subset of M. This shows that we can apply the same technique as in the proof 
of Theorem 4: shift the contour of integration, collect all the residues and estimate the resulting 
integral. The details are left to the reader. □ 

Remark 5. Note that if p = pi G (0, 1) then the coefficient bi^i defined by (4.14) is strictly negative. 
Theorem 6 shows that as a; — 0"'' we have 

k(x) = k(0) + k'(0)x + + o(x^+P), 

1 + p 

which implies that in this case k"{0) does not exist. 

We would also like to point out that if conditions (4.11) are not satisfied, then A^^ o-(s) will have 
multiple poles. This is not a big problem, but it implies that the asymptotic expansions (4.16) will 
contain terms of the form ln(x)^, where —a is the pole of A^^^o-(s) and /c is a non- negative integer 
which is not greater than the multiplicity of the pole —a. Also, results similar to Theorem 6 can be 
derived for a more general class of Levy processes, for example for processes which have jumps of 
rational transform, see [19] for results in the case a = 0. 
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